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^ Abstract 
g 

I The paper concerns fractal homeomorphism between the attractors of two bi- 

afHne iterated function systems. After a general discussion of bi-affine functions, 
conditions are provided under which a bi-affine iterated function system is con- 
CD tractive, thus guaranteeing an attractor. After a general discussion of fractal 
homeomorphism, fractal homeomorphisms are constructed for a specific type of 
bi-afHne iterated function system. 

o 

^ 1 Introduction 

^ ^ The purpose of this paper is to investigate bi-affine iterated function systems and fractal 

^ homeomorphisms between their attractors. The class of bi-affine functions from to 

^ is more general than affine transformations but less general than quadratic transfor- 

mations. These are the functions / : — > that are, for a fixed x or a fixed y, affine 
in the other variable: 

/((I - + axa,!/) = (1 - a) /(xi,?/) + a/(a;2,i/) and 
/(x, (1 - a)yi + ayi, ) = (1 - a) /(x, yx) + a f{x, 7/2) 

for all xi,X2,yi,y2,C( e R. Interpreted geometrically, these equations mean that 

1. Horizontal and vertical lines are taken to lines, and 

2. proportions along horizontal and vertical lines are preserved. 
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This elementary class of functions, with connections to classic geometric results of Bri- 
anchon and Lampert dating back to the 18th century, proves extremely versatile for the 
applications described in this paper. 

Our main motivation for investigating bi-affine functions comes from the representa- 
tion and transformation of certain fractal images. A standard method for constructing 
a deterministic self-referential fractal is by an iterated function system (IFS). The at- 
tractor of the IFS is usually a fractal. Barnsley [21 [3] has introduced a method for 
transforming the attractor of one IFS to the attractor of another IFS, a method that 
has applications to digital imaging such as image encryption, filtering, compression, wa- 
termarking, and various special effects. Figure [T| explained in more detail in Section |5} 
is obtained by applying such transformations - called fractal homeomorphisms. In con- 
structing a fractal homeomorphism it is convenient to use an IFS whose maps have nice 
geometric properties but that are not too complicated. Linear transformations are easy 
to work with and have the property that lines are taken to lines. Affine transformations 
are not much more complicated than linear transformations and do not have the re- 
striction that the origin be taken to the origin. There is a tradeoff; the more properties 
required, the more complicated the mapping. In the fractal geometry literature, in par- 
ticular for fractals constructed using an iterated function system, affine transformations 
are frequently used. For the applications described in Section [5| requiring slightly more 
general functions, bi-affine functions are ideally suited. 




Figure 1: Two fractal homeomorphisms applied to the original picture. 



The paper is organized as follows. The geometry of bi-afiine functions is the subject 
of Section [2j Theorem [T] gives basic properties of a bi-affine function, in particular prop- 
erties of the folding line and folding parabola. Theorem |2] gives a geometric construction 
for finding the image of a given point under a bi-affine function and precisely describes 
the 2-to-l nature of a bi-affine function. Section |3] provides background on iterated func- 
tion systems and their attractors. It is a classic result that, if an IFS is contractive, then 
it has an attractor. Theorem [3] gives fairly general conditions under which a bi-affine 
IFS is contractive. Fractal homeomorphism between the attractors to two IFSs is the 
subject of Section |4j The construction of a fractal homeomorphism depends on finding 
shift invariant sections of the coding maps of the two IFSs. Theorem |4] states that any 
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such shift invariant section comes from a mask. (The terms coding map, section, shift 
invariant, and mask are defined in Section |4j) Theorem |5] concerns how to obtain a 
fractal homeomorphism between two attractors from the respective sections. The visual 
representation of a fractal homeomorphism for bi-affine IFSs is the subject of Section |5] 
Theorem [6] states that a particular type of bi-affine IFS can be used to construct fractal 
homeomorphisms. The main issue is proving the continuity of the map and its inverse. 
The pictures in Figure [T] are obtained by this method. 



2 Geometry of Bi-affine Functions 

Boldface letters represent vectors in M^. A function / : — )■ is called bi-affine if it 
has the form 

/(x, y) = a + hx + cy + dxy. (2) 

It is easy to verify that the class of bi-affine functions is exactly the class characterized 
by the two properties listed in the introduction. The class of bi-affine functions is not 
closed under composition, but the composition of a bi-affine and an affine function is 
bi-affine. 

Call a bi-affine function non- degenerate if d 7^ and neither b nor c is a scalar 
multiple of d. In particular, neither b nor c is the zero vector. If d = 0, then / is 
"degenerate" in the sense that it is affine and well understood. If d 7^ and both b 
and c are scalar multiples of d, then / is "degenerate" in the sense that the range of 
/ degenerates to a line. If d 7^ and just one of b and c is a scalar multiple of d, 
then / is "degenerate" in the sense that the image of the folding line, as defined in the 



next section, is just a point, a fact that can be verified by equation ^ in the proof of 
Theorem [H 

Basic properties of bi-affine functions are described in this section. According to 
statement 2 of Theorem [T] below, the image of a line L under a bi-affine function is a 
parabola. Such a parabola can be degenerate in the sense that it is either a line (focal 
distance 0) or a line that doubles back on itself (focal distance 00). The first case occurs 
if and only if L is parallel to either the x or y-axis. An example of the second case is 
(X, Y) = f{x, y) = {—x + xy, —1 — y + xy), in which case the image of the line y = x + 1 
is given by the parametric equation X = t'^, Y = t"^ — 2, a. degenerate parabola with 
vertex at (0, —2) that doubles back along the line y = x — 2. 

For vectors a and b, let |a b| denote the determinant of the matrix whose columns 
are a and b. For a non-degenerate bi-affine functon /, call the line Lf with equation 

|b d| s + |d c| ?/ = |c b| (3) 

the folding line of the function /. Note that, by the non-degeneracy of /, neither |b d| 
nor |d c| is zero, and hence the folding line is not parallel to either the x or y-axis. The 
terminology "folding line" is justified by the next theorem, in particular statement 4. 
The image Pf = f{Lf) of the folding line is called the folding parabola, which, according 
to statement 3 is non-degenerate. Let and L_ denote the closed half spaces above 
and below Lf, respectively. 
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Theorem 1. If f is a non-degenerate bi-ajjine map, then 



1. any line parallel to either coordinate axis is mapped to a line; 

2. any line is mapped to a, possibly degenerate, parabola; 

3. the folding parabola Pf := f{Lf) is non- degenerate; 

4- the map f is injective when restricted to either or L_. 

Proof. Statement (1) follows from the fact that, for a fixed a; or a fixed y, a bi-afiine 
function is an afiine function in the other variable. 

The image of a line under a bi-affine transformation has a parametric equation of 
the form f{t) = u + vt + wt^, which is a, possibly degenerate, parabola. 

To show that Pf is non-degenerate, first note that the image of the line Lf under / 
is given by the parametric equation (X(x), Y{x)) with parameter x by 

(X,Y) =f(x,x) = — ^ [(Idcl a+ Icbl c) + (Idcl b - Ibdl c + Icbl d)x + Ibdl dx^l 
|dc| 

= p'— I [(|dc| a+ |cb| c) + 2|db| ex + |bd| dx^] , 

(4) 

the first equality simply by substituting from equation 3_ into equation [2] and the second 
equality by a direct calculation. The tangent vector to this parabola is 

^/ ^ 2|dbL , , 
T(x) = ,' ,' c -dx). 
|d c| 

By the non-degeneracy of /, the vectors c and d are linearly independent, implying 
that the the direction of the tangent vector is not constant. Hence the parabola is not 
degenerate. 

The Jacobian determinant of a bi-affine function is |bd|x + |dc|?/ — |cb|, which is 
nonzero except on the folding line. It follows from the inverse function theorem that / 
is injective when restricted to either L+ or L_. □ 

A parabola P divides the plane into two regions; let P denote the closed region 
"outside" P (and including P). The set Pf will be called simply the parabolic region of 
f. For a point p = {x,y) G M^, let 



|cb 


— 1 


d c 


y 




b d 





|c b| — |b d| X 
\dc\ 



(5) 



This somewhat complicated formula is merely an analytic expression of the simple ge- 
ometry shown in Figure |2j Statement 3 in Theorem |2] below makes precise the 2-to-l 
nature of a bi-affine function given in statement 4 of Theorem [Tj Statement 1 in Theo- 
rem |2] is illustrated in Figure |3| Statement 2 gives a geometric construction of the image 
of a given point under a bi-affine function and is illustrated in Figure |4j 
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Figure 2: Folding Line Lf : /(p) = p* . 




Figure 3: Folding parabola P/, parabolic region Py, tangent line Ta- 




Figure 4: The image of a point p under a bi-affine function. The point /(p) has 
"coordinates" {A,B). 

Theorem 2. Assume that f is a non- degenerate hi-affine map. We use the notation La 
for the line x = a and Lf, for the line y = b. 

1. If A = {a,b) is any point on the folding line Lf and Ta is the tangent line to Pf 
atfiA), thenTA = f{La) = f{L,). 

2. Ifp = {a,b) e ]R^ then /(p) = TaHTb, where A = LfHLa and B = Lf n Lb. 

3. f{L_^) = f{L_) = Pf and, in particular, /(p) = /(p*) for all p = (x, y) G M^. 
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Proof. Concerning statement 1, any horizontal or vertical line L intersects in a single 
point. Therefore f{L) intersects f{Lf) = Pj- in a single point, which implies that f{L) 
is tangent to Pf. Since the intersections Pj fl f{La) and Pf fl f{Lb) both consist of the 
same single point, f{La) and /(I/ft) both equal the tangent line T4 to Pf at f{A). 

By statement 1, the point /(p) lies both on the tangent to Pf at f{A) and on the 
tangent to Pf at f{B). This proves statement 2. 

Consider p and p* as in Figure [2] Statement 2 implies that /(p) = T^nTg = /(p*). 
In particular f{L^) = f{L^). Since the union of all tangents to Pf is the parabolic 
region Pf of /, we have f{L+) = f{L^) = /(M^) = p^. □ 

It is a consequence of Theorem [2] that the parabolic region can be coordinatized as 
follows. Each point p G P/ has a unique set {A, B} of (unordered) coordinates where 
A and B are points on the folding line. Specifically p = fl Te. This is illustrated in 
Figures |4] and [5] 




Figure 5: Coordinatization of Pf. Note that A- = f{Ai), 1 < i < 11, and p' = /(p). 
The point p' has coordinates {A3, A9}. 

The folding parabola Pf can be constructed geometrically as follows (See Figure [6]). 
Choose four pair wise distinct points Ai, A2, A3, A4 on the folding line and construct the 
four tangent lines T^^ , 5 5 • It is a direct consequence of the classic Brianchon 
Theorem (actually the converse) that there is a unique parabola with these lines as 
tangents. According to what has been shown above, this must be the folding parabola. 
The parabola can be explicitly constructed using a theorem of Lambert. According 
to Lambert, the circumcircle of a tangent triangle of a parabola (see Figure 6) goes 
through the focus of the parabola. So the focus F is determined as the intersection of 
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three such circumcircles. Then reflect F about two of the tangents to get two points on 
the directrix. 




Figure 6: Construction of the folding parabola: focus F and directrix D. 



3 Iterated Function Systems 

This section reviews the standard notation and definitions related to iterated function 
systems (IFS). These concepts are then applied to the bi-affine case. 

Let X be a complete metric space. If : X — )■ X, m = 1, 2, . . . , M, are continuous 
mappings, then J-" = (X; /i, /2, /a/) is called an iterated function system (IFS). An 
iterated function system that consists of bi-affine functions will be called a bi-affine IFS. 
To define the attractor of an IFS, first define 

HB) = U f{B) 

for any 5 C X. By slight abuse of terminology we use the same symbol J-" for the IFS, 
the set of functions in the IFS, and for the above mapping. For i? C X, let J^^{B) 
denote the /c-fold composition of J-", the union of /j^ o /^^ o ■ ■ ■ o fi^{B) over all finite 
words iii2 ■ ■ - ik oi length k. Define J^{B) = B. A nonempty compact set A C X is said 
to be an attractor of the IFS J-" if 

1. T{A) = A and 

2. limk^oo J^^{B) = A, for all compact sets 5 C X, where the limit is with respect 
to the Hausdorff metric. 

Attractors for bi-affine IFSs consisting, respectively of 2, 3 and 4, functions are shown 
in Figure [7| 
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Figure 7: Attractors of IFSs consisting of two, three, and four bi-affine functions, re- 
spectively. 

A function / : X — > X is called a contraction with respect to a metric d if there is an 
< s < 1 such that d{f{x), f{y)) < s d{x, y) for all x,y E M". An IPS with the property 
that each function is a contraction will be called a contractive IFS. In his seminal paper 
Hutchinson ^ proved that a contractive IFS on a complete metric space has a unique 
attractor. Theorem [3] below gives fairly general conditions under which a bi-affine IFS 
is contractive. 

Let □ denote the unit square with vertices (0,0), (1,0), (1, 1), (0, 1). The shape of 
the image of □ under a non-degenerate bi-affine map / depends on the location of the 
folding line Lf relative to □. It follows from Theorem |2] that there are three possible 
cases as shown in Figure [s] It is only in Case ^ {Lf disjoint from the interior of □) that 
the image of the sides of □ form a convex quadrilateral. If this is the case, call / proper. 

Note that the unique bi-affine function taking (0, 0), (1, 0), (1, 1), (0, 1) to the points 
Po,Pi,P2,P3, respectively, is 

fix, y) = Po + (pi - Po)a; + (ps - Po)y + (P2 + Po - Pi - P3)xy. (6) 

Theorem 3. Let f{x, y) = Po + (pi - Po)x + (ps - Po)y + (P2 + Po - Pi - P3)xy be a 
proper, non- degenerate bi-affine function. If there is an s,0 < s <1, such that 

\Pi+i-Pi\ < s, \Pi+2-Pi\ <V2s, \pi+i + pi^i - 2pi\ < V2s 

for i = 0, 1, 2, 3 [mod 4), then f is a contraction on □. 

In terms of the quadrilateral P0P1P2P3, the first of the three inequalities states that 
each side has length less than or equal to s, the second that each diagonal has length 
less than or equal to \/2s, and the third that vector sum of any two incident sides has 
length less or equal to \/2s. For a bi-affine function taking □ into itself, for example, 
these conditions are not too restrictive. Two lemmas help in proving Theorem [3) 

Lemma 1. If f is a proper, non- degenerate bi-affine function, then f is injective when 
restricted to □. 
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Case 1 





Case 2 



Case 3 




Figure 8: The image of the sides of the unit square (on the left) is shown by thick hues 
(on the right): A = /(A), B' = /(S), C = f{C), D' = f{D). 



Proof. By the comments above, the folding line L f lies outside the interior of the square 
□. The lemma then follows by statement 4 of Theorem [1} □ 

Lemma 2. Let f be a non- degenerate hi-affine function that is injective on and let 
W O \z\ be any rectangle with sides parallel to the x and y-axes and with diagonal of 
length piW). Further let piiW) and p2{W) be the lengths of the two diagonals of f{W) 
and let 

M{po) = max max{pi(W) , p2iW)) . 

W:p{W)=po 

For any po < \/2, if Wq is a rectangle that maximizes M(pq), then Wq and □ have a 
common vertex. 

Proof. If one side of Wq lies on the line x = or x = 1 and another side lies on y = 
or y = 1, then the proof is complete. So, without loss of generality, assume that Wq 
has no side that lies on x = or a; = 1. Let Vq be the rectangle bounded by the 
lines x = 0,x = 1 and the lines determined by the upper and lower sides of Wq. Let 
A, B, C, D be the vertices of /(Vq). Then by the conditions 1 in the Introduction, there 
is an a such that the four vertices of f{Wo) are (1 — a)A + aB, {1 — a — A)A + (a + 
A)B, (1 - a - A)C + {a + A)D, (1 - a)D + aC, where A is the horizontal length of Wo- 
As a varies between and ao, the rectangle Wq shifts left or right, from the extreme 
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left side of □ to the extreme right side of □. The lengths of the two diagonals of /(Wq) 
are \Ui + (B + C - A - D)a\ and \U2 + {B + C - A - D)a\, where vectors Ui and 
U2 depend on A,B,C,D and A. As a varies in the range < a < ao, the quantities 
Ui + {B + C — A — D)a and U2 + {B + C — A — D)a describe (parallel) line segments. 
Hence the maximum of \Ui + {B + C — A — D)a\ and \U2 + (5 + C — A — D)a\ occur 
at an end, i.e. a = or a = ao, contradicting the assumption that Wq has no side that 
lies ona; = Oora; = l. □ 



Proof, (of Theorem |3]) By the continuity of / and the compactness of □, it is sufficient 
to show that there is an s, < s < 1 such that \ f{x,y) — f{x',y')\ < s \{x,y) — (x',y')|. 
Let Wq be the rectangle whose diagonal is the line segment joining (x, y) and {x'y'). By 
Lemma [l| the function / is injective on □, and by Lemma [2| the quantity 

D \f{x,y) - f{x',y')\ w / / /A 
\{x,y) - {x',y'}\ 

is maximized when Wq and □ have a common vertex, i.e. when Wo lies on the corner of 
□. Without loss of generality, it may be assumed to be the lower left corner. Otherwise, 
replace / with the composition of / with the rotation that moves the lower left corner 
to the relevant corner. Now let Ax = \x' — x\ and Ay = \y' — y\. With b = pi — po, c = 
Ps ~ Po; d = P2 + Po — Pi — Ps and setting r = we have two possible formulas for R: 



,2 |/(Ax,At,)-/(0, 



or 

|/(Ax,0)-/(0,Ay)p _ 1 

Ax,A?/)|2 l + r2' ' ■ 

Since < Ax < 1, the quantity R^ in the first formula is maximized when either Ax = 
or Ax = 1. So it is now sufficient to prove that 

max ( -lb + rcP, -lb -|- rc -|- rdP, -|b — rcP) < 1. 

The inequalities in the statement of the theorem, in terms of b, c, d are that |b|, |c|, |b-|- 
d|, |c -|- d| are all less than or equal to s and |b — c|, |b -|- c -|- d|, |b -|- c|, |b — c -|- d|, |c — 
b + d|,|b-|-c-(-2d| are all less than or equal to a/2s. Maximizing the three quantities 
over r, in order to verify the above inequality, is a not quite trivial calculus problem 
whose details are omitted. □ 

Example 1. Consider a bi-affine IFS J-" = {□; /i, /2, /s, /4}, where the four functions 
are determined by the images of the four vertices of □ as shown in Figure |9} Each of 
the four images /»(□), i = 1, 2, 3, 4, of the square contains exactly one vertex of □. The 
attractor of J-" is □ itself. For "most" choices of the center and side points, J-" satisfies 
the conditions of Theorem [3] and hence J-" is a contractive IFS. Why this simple example 
should be of interest is the subject of the next two sections. 
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Figure 9: An IFS consists of four bi-affine functions. Points labeled with lower case 
letters are the images of points labeled with upper case letters. The attractor is the unit 
square. 



4 Fractal Homeomorphism 

In the case of a contractive IFS, it is possible to assign to each point of the attractor 
an "address". Given two IFSs J-" and Q with respective attractors Ap and Aq, a fractal 
homeomorphism is basically a homeomorphism h : Ap ^ Aq that sends a point in 
to the point in Aq with the same address. To make this notion precise, let J-" = 
(X; /i, /2, /at) be a contractive IFSs on a complete metric space X with attractor A. 
Let VL = {1, 2, . . . , N}°° denote the set of infinite strings using symbols 1, 2, . . . , A^, and 
for a & Q, let cr\k denote the string consisting of the first k symbols in a. Moreover, if 
a = toil, "^2 ■ ■ ■ , then we use the notation 

fa\„ ■= /io ° /ii ° ■ ■ ■ ° /i„- 
Now define a map n : Q ^ A, called the coding map, by 

7r(a) := lim 

fe— >oo 

It is well known [H H] that the limit above exists and is independent of x G X. Moreover 
TT is continuous, onto, and satisfies the following commuting diagram for each n = 
1,2,. ..,iV. 



n 




n 


71 I 




i 71 


X 




X 




fn 





The symbol s„ : — ?■ 1) denotes the inverse shift map defined by s„(cr) = na. A section 
of the coding map vr is a function r : Q ^ A such that tt o r is the identity. A section 
selects, for each x G A, a single address in Q from the ones that come from the coding 
map. Call the set fir '■= '^{^) the address space of the section r. Let S denote the 
shift operator on fl, i.e, S{na) = a for any n G {1, 2, . . . , A^} and any a E fl. A subset 
W G fl will be called shift invariant if a E W implies that S{a) G W. If fir is shift 
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invariant, then r is called a shift invariant section. The following example demonstrates 
the naturalness of shift invariant sections. 

Example 2. Consider the IFS J-" = (M ; /o, /i) where /o(a;) = and /i(a;) = | x + |. 
The attractor is the interval [0, 1]. An address of a point x is a binary representation of 
X. In choosing a section r one must decide, for example, whether to take r(|) = .01 or 
r(|) = .00111 ■■■ . If the section r is shift invariant, this would imply, for example, that 
if r(|) = .00111 ■ ■ ■ , then r(i) = .0111 ■ ■ ■ , not r(^) = .100 • • • . 

Lemma 3. With notation as above, a section r of an IFS is shift invariant if and only 
if, for any x E A, if t(x)|i = n, then {S o r)(x) = (r o 

Proof. Given the right hand statement above, we will prove that r is shift invariant. 
Assume that na G Qr- Then there is an a; G A such that r(x) = na, and hence 
a = {So t){x) = (r o = T{f-'){x)). Thus a G fi.. 

Conversely, assume that r is shift invariant. Assume that x E A and r(x) = na for 
some 0" G fi. By shift invariance, there is a y G A such that T{y) = a. Now 

a; = (tt o r)(a;) = 7r(no-) = /„( lim /<,|,(A)) = /„(7r(r(y))) = /„(?/). 

fe— >-oo 

Therefore ?/ = fn^{x) and 

(5 o r)(a;) = S{na) = a = T{y) = T{r~\x)) = (r o f~^){x). 

□ 

Call an IFS injective if each function in the IFS is injective. Theorem |4] below states 
that every shift invariant section of an injective IFS can be obtained from a mask. For 
an IFS J-' with attractor A, a mask is a partition M = {Mi, 1 < i < N} of A such 
that Mi C fi{A) for all /j G J-". Given an injective IFS J-" and a mask M, consider the 
function T : A ^ A defined by T(x) := when a; G Mj. The itinerary tm{x) of a 

point a; G v4 is the string io ^2 ■ ■ ■ G fi, where ik is the unique integer 1 < < such 
that 

T\x) G M,,. 

Theorem 4. Lei be a contractive and injective IFS. 

1. If M is a mask, then tm is a shift invariant section ofn. 

2. If T is a shift invariant section of tt, then r = T]\j for some mask M. 

Proof. To show that tm is a section, let x E A. Then (tt otm){x) = limfc_>oo fTM(x)\ki^)- 
It follows immediately from the definition of tm that x G /rM(2:)u(^jfc+i) — frM{x)\ki^) 
for all k. Hence (tt o tm){x) = x. Concerning the shift invariance, it follows from the 
definition of T that the following diagram commutes. 



A 


T 
-> 


A 








n 


S 


n 
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If 0" G f^TM5 then there is an x G A such that a = Tm{x) and, from the diagram, 
S{a) = (S oTM)ix) = {tmoT){x) = TuUn^ix)) G Vt^j^ for some n G {1,2, . . .,N}. 
Concerning the second statement, define a mask M = {Mi, 1 < i < N} as follows: 

Mi = {x : t{x) = 1(7 for some a G Q}. 

It is sufficient to show that Mj C fi{A), and that tm{x) = t{x) for all x E A. If a; G Mj, 
then r(x) = ia for some a G and a; = (vr o r)(x) = /^(limfc^oo /o-|fe(^)) ^ 

To show that rM(x) = t{x), let r(x) = joji j2 ■ ■ ■ and rM(a;) = kokik2---. That 
Jo = ^0 follows form the definitions. By induction, assume that ji = kii = 0,1, ... , m—1. 
Applying Lemma [3] for m times yields 

Jm = o r){x))U = r{f.^_, o . . . o 41 o f7^\x)\„ 

where the jVs are determined by the recursive formula T{fJ^\^ o ■ ■ ■ o fr^ o fjQ^{x)\i = jr- 
By the definition of the mask, T{f~^_^ o ■ ■ ■ o o /^^(x)) G Mj^^. But by the definition 
of the itinerary, the /cj's are determined by the recursive formula 

r{U.\.o...of-^of-^{x))eM,,. 

ioT r = 0,1,2, ... ,m. But, since ki = jii = 0,1, ... ,m — 1, we have f{fj^_^ o ■ ■ ■ o fr^ o 
/jQ^(x)) G Mfc^. Therefore fc^ = j^. □ 

To define fractal homeomorphism, consider two contractive IFSs J-" = (X; /i, /2, /at) 
and Q = {li;gi,g2, ...^Qn) with the same number N of functions on a complete metric 
space X. Let Ap and Aq be the attractors and t^f and ttg the coding maps of J-" and 
^, respectively. A homeomorphism h : Ap ^ Aq is called a fractal homeomorphism if 
there exist shift invariant sections Tp and such that the following diagram commutes: 

Ap Ag 
rF\ ' /TG (9) 

i.e., the homeomorphism h takes each point x G with address o = Tp{x) to the 
point y G Ag with the same address a = TG{y). Theorem |5] below states that the 
fractal homeomorphisms between attractors Ap and Ag are exactly mappings of the 
form ttg o Tp OT 7Tp o tg for some shift invariant sections Tp, tg. 

Theorem 5. Let T and Q he contractive IFSs. With notation as above: 

1. If h : Ap — 7- Ag is a fractal homeomorphism with corresponding sections Tp and 
Tg, then Qrp = ^tq- Moreover h = ttg ° Tp and h^^ = irp o tg ■ 

2. If Tp is a shift invariant section for and h := ttg^tp is a homeomorphism, then 
h is a fractal homeomorphism. 
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Proof. Concerning statement 1, since /i is a bijection, the commuting diagram [9^ implies 
that the images of Tp and tq are equal, i.e., Qrp = ^tg- Now Tp = tq o h from the 
diagram implies ttg o Tp = ttg ° Tp = [ttg o tg) o h = h. The formula involving is 
likewise proved. 

Concerning statement 2, the section Tp is a bijection from Ap onto ^rp- Since /i is 
also a bijection, the equality h = ttg o Tp implies that tcgIuf^ the restriction of ttg to 
is a bijection onto A^. If tg is the inverse of ttgIq^, then tjt' and tg satisfy the 



commuting diagram 9 , That Tp is shift invariant means that Qt-q = , i-e., tg is 



shift invariant. □ 



5 Image from a Fractal Homeomorphism 

This section concerns images on the unit square □. Define an image as a function 
c : □ — 7- C, where C denotes the color palate, for example C = {0, 1, 2, . . . , 255}'^. If h is 
any homeomorphism from □ onto □, define the transformed image h{c) : □ — > C by 

h{c) := CO h. 

We are interested in the case where /i is a fractal homeomorphism. The remainder of 
this section concerns fractal homeomorphism based on bi-affine IFSs with four functions 
as described in Example [1} 

Consider Example [T] depicted in Figure |9] For the bi-affine IFS 
J-" = {□; /i, /2, /a, f^}, we will construct a section Tp that is referred to in |2] as the top 
section. Consider the mask Mp = {Mi, M2, M3, M4} defined recursively by 

M, = /,(n)\U/,(n) 

for i = 1,2,3,4. Explicitly, Mi is the closed quadrilateral ATOQ, M2 is the open 
quadrilateral OQBR together with the segments {Q,B], [B,R], [R,0), M3 is the open 
quadrilateral ORCS together with the segments {R,C],[C, S],[S,0), and M4 is the 
open quadrilateral OSDT together with the segments [S^D], [D^T). The section Tp 
corresponding to the mask Mp is given by t{x) = max7r~^(x), where the maximum is 
with respect to the lexicographic order on Q. 

Now consider a second bi-affine IFS Q = {O; gi, g2, ds, 94} of the same type with 
points O' ,Q' , R' , S' ,T' replacing 0,Q, R, S,T, and with mask Mg defined exactly as 
it was for Mp. The masks Mp and Mg induce shift invariant sections Tp and tg, 
respectively, as verified by Theorem |4j Theorem |6] below states that ttg o Tp and upotg 
are continuous and hence, by Theorem |5| fractal homeomorphisms. 

To prove Theorem [6} the following lemma will be used. The proof is routine and will 
be omitted. All partitions V will be of the unit square □ into regions whose closures are 
topological polygons. The dual graph of such a partition is the graph T-p whose points 
are the regions and where two vertices are joined if and only if the corresponding regions 
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share a side. A partition V is nested in partition Q if each region in V is contained in 
some region of Q. Assume that partition Vi is nested in partition V2 and Qi is nested 
in Q2, and that there are graph isomorphisms $1 : V-p^ — ^ Tg^ and $2 : Tp^ Vq^. 
Call $1 and $2 compatible if whenever Pi G Vi and P2 G V2 with Pi C P2 we have 
$i(Pi) C ^2{P2)- The mes/i l^l of a partition P is the maximum diameter of the 
regions. If lim„^oo I'^nl = 0, then, for any s G □, there is a unique nested sequence {Pn} 
of regions P„ G Vn such that x = HneN -^n- 

Lemma 4. Lei Vn and Qn, n = 0,1,2, ... , be two nested sequences of partitions of the 
unit square □ with lim„_5.oo |'Pn| = lini„_s>oo \ Qn\ =0. Assume that there are corresponding 
sequences of compatible graph isomorphisms : G-p^ — >■ Gq^ . The map /i : □ — □ 
defined as follows is a homeomorphism. For a; G let x = HneN -^n £ □ ^^/^ -Pn £ 
and define h{x) := f]^^j^^niPn) ■ 



Figure 10: The dual graph of the partition V2- 



Theorem 6. For the two bi-affine IFS's T = {□; /i,/2,/3,/4} and 
Q = {D] gi, 92, 93, 94} defined above, the map h = -kq ° tf is a homeomorphism. 

Proof. For each n > 0, let f2„ denote the set of strings of length n using symbols 
{1,2,3,4} For the IFS J-', define a partition Vp = {Pa '■ cr G fi„} of □ recursively by 
taking Vp = Vf and 

= {Po, =P.^ fAP,) : (renn,l<J< 4}. 
A straightforward induction shows that {Mp} is a nested sequence of partitions of 



□. The dual graph F'^ of Vp is the grid graph shown in Figure 10 for n = 2. This 
construction of a nested partition can be repeated for the IFS Q. Since the obvious 
graph isomorphisms between F^ and Tq are compatible with the nested partitions. 
Lemma |4] implies that h is a homeomorphism. □ 
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